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The effects of boundary perturbations on finite-amplitude 
acoustical standing waves in rectangular, rigid-walled 
Cavities were investigated using non-linear theory. When a 
high amplitude standing wave of frequency w is generated in 
a cavity, non-linear effects will cause a stimulation of 
certain normal modes whose resonance frequencies are integer 
multiples of w. Previous experimental observations revealed 
that there could be excitation of other normal modes, not 
belonging to the family of the driven mode, which was not 
predicted by the non-linear theory. 

The purpose of this research was to investigate the 
possibility that deviations from the idealized geometry could 
account for these observations. Of the various mechanisms 
possible, this work investigated the possibility of these 
unpredicted excitations occurring through a non-linear 
mechanism, 

The standing waves that ex1st 1n an ideal cavity must be 
corrected when the boundaries are irregular. The non-linear 
interaction between these standing waves and the corrections 
was studied. The ability of this interaction to excite 
standing waves other than those predicted in the ideal case 
was verified. A specific example was worked out demonstrating 
an unpredicted excitation, the strength of which was on the 


order of the magnitude of the boundary perturbation parameter. 
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ieee ORUCT ION 


The purpose of this research was to investigate some 
of the effects of boundary wall perturbations on finite- 
amplitude standing waves in a rectangular cavity. The 
investigation was prompted by an examination of the experi- 
mental results of Coppens and Sanders [1] and the research 
of DeVall [2], which suggested the existence of the 
excitation of modes other than those belonging to the family 


of the driven mode. 





PE bacnwCROUND 


An acoustical wave traveling in a medium will have its 
form distorted as predicted by the non-linear hydrodynamic 
equations. If the medium is absorptive then only waves of 
relatively high amplitude will develop distortion before 
becoming completely attenuated due to the absorptive process. 

At the United States Naval Postgraduate School, Coppens 
and Sanders [1], DeVall [2], and others [3,4] have dealt 
with the study of finite amplitude waves in rigid-walled 
rectangular cavities and tubes. One interesting result of 
these cavity experiments was the appearance of excitations 
Sremedes witch were not family members of the driven mode. 
For example, assume that the cavity is driven in the (0,1,0;1) 
mode, which indicates the excitation is of the form 


cos KY cos wt , where 


and 


Then a non-familial mode is any which does not have the 
fourier oOs;n) . Lt was proposed [5] that the preseneec 


of these non-familial modes could be attributed to the cavity 


AG 





walls being non-ideal (not perfectly flat) and that perhaps 
the previous small disparity between theoretical predictions 
and experimental observations might be partly explained by 
the presence of these non-familial modes. 

The non-linear wave equation applicable to this problem 


can be expressed [6] as 











2 U 
24 P ee 24 all 2 Ge? 
Sor Er. . 250 aaa Aone 
Pog ot hed 0 
(22s) 
U 
ie. Same Dea 0.2 
ts Cy «(UV ee a Gee 
PoXo 


oO is defined as 


where L represents the loss terms appropriate to each 
frequency component in P . The RHS of the non-linear 
equation may be interpreted as a forcing function (consisting 
of both classical and non-linearly generated terms) which, 

in first-order perturbation theory, is the classical solution 
for the pressure. The first-order perturbation solution may 
be considered to describe the effects of self~interaction of 
the classical solution, whereas higher-order perturbation 


solutions describe the effects of interaction of the 
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Seeeme mon-linear Solution with itself. If the cavity is 
driven at frequency w, non-linear terms force the 
existence of all terms containing the form nw . Each such 
term whose frequency is close to the resonance frequency of 
a standing wave of the cavity and whose spatial behavior 
matches that standing wave can strongly excite that standing 
wave to a degree determined by QQ, the guality factor of 
the resonance, and by the difference between nw and the 
resonance frequency of the standing wave. In a rectangular 
cavity, it can be seen [1] that if the classical solution 

is the (msm, ,m,7q) mode, then the non-linear wave equation 
will force the strong excitation of all modes 

(nm, nm, -nm,;nq) which constitute the family. Any mode of 
frequency nw which is not a member of the family should 
not be strongly excited. The purpose of this research is to 
see if the presence of wall irregularites can explain how 


non-family members may be strongly excited, and to present 


an example to support this theory. 


12 





Tie. THESNON-LINEAR TP ROL EE 


Fae JSCAVITY DESCRIPTION 


Assume a perfectly rigid-walled rectangular cavity has 
one wall perturbed such that the cavity dimensions are 


Lf1+sf(y) 1, L L. A top view depicting the XY plane 


eee 2 
and showing the boundary irregularity 1s given in Fig. 1 


below. Also assume the perturbation on the boundary 1s 


small compared to the cavity dimensions, |éf(y)| << 1 
y 
L 
ys 
ig x 
x 
FIGURE 1 


The cavity is to be excited by a source near one of the 


corners in such a way that the (n,m,%;q) mode is excited. 


is 





B. CLASSICAL SOLUTION FOR IDEAL BOUNDARIES 


For a rigid-walled rectangular cavity with ideal 


boundaries (6=0), the classical linear solution for the 


pressure Pp US subject tO une, LOoltowing Condi cions, 


namely 
2 
[TPp = 9 
and VP9 "2 = 0 aie x=0,L, 
ae 
z=0,L, - 
where n is the local normal to the ideal surface. 


0 


solution for Po in terms of Mach number is 





Pp 
= M cos kx cos KY cos k,2 cos (wu t+9) 

210 
where 

=. aq 

ky, = L, 

m7 

kK = 

y Ly 

= i 

ky L, 
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i -Srsade) 


3.2) 


The 


(335) 


(3.4) 





and 


(Ay ¢ = ko = ve 44k 44h 7 (35) 
VY Z 


where we is the frequency associated with the (n,m,2%;q) 


mode. 


Geo tAoo CAL SOLUTION FOR REAL BOUNDARIES 

Since the cavity is assumed to have one perturbed wall, 
1t is reasonable to expect the solution for reai boundaries 
to differ from that for ideal boundaries. The classical 
linear solution for real boundaries, Py , may then be 
considered a summation of the classical linear solution for 


ideal boundaries plus perturbation correction terms due to 


the irregular boundary: 


= 2 


As long as the magnitude of the boundary perturbation is kept 
small, terms higher than first-order can be considered 


imscgmratcant, so that 
Po = Po + Op, (EO Tarse Order) (3.6) 


ana must be a solution to 


Cl Po = QO (3% 7) 


E5 





with boundary conditions 


ean = 0 at x=0, L[1+6f(y)] (3.8) 
y=-0, oe 
z=0, L, ’ 


“Aw 


where n, the local normal to the real surface, is obtained 


by taking the divergence of the equation for the perturbed 





boundary, 1.e.,; es] Uae (yn) . Thus, to first order 
ieee Oy 

nan we “A - gt (yy) “A 

n= x L, 6 By a 


and when the above is used in (3.8) the result is 


Py — 2 8 Po Jf (y) = 0 
o* Mxeb(lest(y)]  * 8Y 8Y [xen pitse(y) ] 





(37579) 


D. THE APPLICABLE NON-LINEAR EQUATION 

A real wave in a real cavity will of course have non- 
linearities introduced as it propagates through the cavity. 
Let Py be the correction term to the pressure due to the 


non-linearities of the medium. The total pressure P is then 


(3200) 


16 





When (3.10) is inserted in (2.1), the LHS of the non-linear 


equation becomes 


jo ag Dy ae 


Cc ae (es 
oe cits ar. a ullie a 








mae Lirst term above is known to be zero from (3.7), so the 
LHS of (2.1) becomes simply 
al 


2 
Se alle ay 
0-0 





Now if standing waves are to be produced in the cavity, Py 
must consist of terms of the form cos kx cos arred cos k,2 cos wt 
Using the basic definitions 


> 
Tn 2 0 
(On Ese 


and 


Vo 


C+ 
H 


Cr+ 


reguires that be cf the form san kx sin ey Sin k 2 Sin wt 
When this information is used, the RHS of (2.1) will consist 
of four terms: 


2 


3 
We moe (cos kx cos KY cos k2 cos wt) 
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_ a : . , 2 
We ae (sin kx Sin KY sin k 2 Sin wt) 


oe ac (Gos kx cos k y cos k.z cos wt) ¢ 
x y Z, 


2 


: . Z 
fee V2 6 «6«CSan kx Sin KY Sin k 2 Sin wt) 


Although the expansion of these will produce many terms, it 
must be noted that only some will produce resonant standing 
waves. All the other terms may therefore be considered to 
give a small contribution when compared to the resonant terms, 
and thus may be neglected. This fact, when used with the 


trigonometric identities 


cosa = > [1 +cos 2q] 
and 
stim” == [1 -cos 2q] 


meauees the RHS of (2.1) to a single term, and after a good 
deal of manipulation the non-linear equation may be 


approximated by 


p , p 
c. 242 sae ere eae Eg (3.11) 
o LIL 2 b. 2 

Poo Poo 








The governing non-linear equation is now in a relatively 
simple form, and Appendix I develops the general form of 


the solution for Py 
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iVe “LES Ee RT URED seOUnDARY 


A. DEVELOPMENT OF A BOUNDARY CONDITION ON Py 

Recall that Py is the first-order correction term due 
to one boundary being irregular. Although it is known that 
for a perfectly rigid-walled cavity the pressure must be a 
maximum at the walls, this fact may not be conveniently 
used with one wall having an irregular form such as 
Gli + of (y)] - it would be highly desirable, therefore, 
to develop a boundary condition on P, to be applied at L, 
This may be accomplished as follows: 

Using the incremental form for a Taylor series expansion 


on Po evaluated at the real boundary ep ae ee Sih) produces 








38P fl a°P 2 
Po| 1 3x L,of (y) + oy 2) IL, 5f (y) } + 
ox = 
x=L xXx=L x= 
x x x 
Sulpstitucion of 


on the RHS, and taking the partial derivative with respect to 


x on both sides gives 


iby, 








ed (mene (lo) 


Following exactly the same procedure for the y coordinate 








gives 
a ane oy De 
y=L, (1+df (y) ] x=L, ie 
2 
,) Po 
+ 5 | le ee 2). ar (4.2) 
oY X=L = 
x 
Recall that (3.9) required 
Ie 
al Ee enon e(y) - 
OX ee A oy 


=L, [1+6£ (y) ] x=L [1+6f (y) ] 


When (4.1) and (4.2) are used in (3.9) the resuit is 


20 
































dp 6 dP d Pp 
—— 5 4 + 0 L 6f£(y) ar 
x=L es x=L Spe x= 
Xx Xx Xx 
Pp 6 dp 3“ 
if 
L be i —: n _ Ts o6ee yale 
x=L Y x=Ly oY x= 
Collecting like-order terms in § then gives 
—— + —— tek (vy) = ah =) 
es x= ox x=L * ee oY x=L 
Xx x 
ore as a boundary condition on Py applied at Lys 
dp dp dP a! 
aX 0 df (y) 0 
aes = L [-fly) —- + ~—- Vie (4.3) 
aes x=L = ox x=L oY Y |x=L 


Pee A FOURLER SERIES REPRESENTATION 

It will be desirable for further development to represent 
the RHS of (4.3) as a Fourier series in cosines, so that Py 
may be expressed aS a Summation of normal modes. First, 


rewrite (4.3) as 


2 
Pp 0 Pg dP 
al 1 = ML [-f fe oily) 70 Aune 
gx J ell (y) wee ay OY <f } (Are 


P2AL 





Bee the RHS of (4.4) equal ‘Scos k 2 cos wt and express it 


as an infinite Fourier series as follows: 


ja co 
S(y) cos Kz cos wt = [ +- LAY cos =) Cos K2Z0COs scr 
L Z, 
N=1 4 
where the An are given by 
My 
2 
AL = ‘6 S (y) cos = ay, N= NOT alsa gecnonene 
y 0 y 


Then (4.4) becomes 


Ay cos Iv cos k_z cos wt. 
L Zz 
If y 


=| 
ee 
| 2 
I 
we 
4- 
ll cr 8 


e N 


mae SOlution for Py is 





Py cos k,'x N 
a ee ee = cos k,z cos wt (EES), 
NG eee Sante 
Py <o x >, Se? 2S Y 


where for N = 0 the coefficient is Ay/2, and Soe 1s such that 


Note that hs is not in the form of a standing wave in the 


cavity required by 


OG 





ee ae (3.4) 


Thus the terms involving Keg in (4.5) must also be expressed 


as a Fourier series: 








B 
“3 + ee a elo Max ; (4.6) 


where the By are given by 


1 


L ' 
5 x cos ky Dt Meee 
Bu ns os if ——- COS TT «Xx, Me="07172, 
S010 ke {Sine xX 
X X "xX 


When the results of (4.6) are inserted into (4.5), Py is 


expressed as a double infinite series: 


= Max Nay 
eo a» Bu cos 7 Ay cos 7 «COS KZ cos wt, (4.7) 
Cc N M x Y 





where again each of the zero coefficients (Aye BD) are to be 
divided by a factor of two. 

The formalism has now been developed whereby given any 
driving mode and boundary irregularity, Py 7 the first-order 
perturbation correction term ~ may be expressed in terms of 


the spatial components of the normal modes of the ideal 


Cava ty. 


Ze 





Ve A SPECLPIC EPAAMPLE 


A. THE COSINUSOIDALLY PERTURBED BOUNDARY 


AS a Specific example utilizing the theory developed 


in the preceeding sections, assume the cavity is driven in 


meme  (0,1,0;1) ~modew “Thus 


Po _ 


a M cos KY cos wt : C5.) 





POSo 


Further assume that the boundary is perturbed cosinusoidally 


such that 


f(y) = cos --y : (iecy 


ime boundary condition on Py then becomes 





ae =) = M - T 7 sin _ sin KY COS mits 
00 x=L Y 


het 


which, when used with the requirement that k, = — ; 
Y 
becomes 
L L 
rs) 1 Se. ull _ Sc odl, 2Ty 
ae Taper) = M Teo k cos wt M = Kk, cos Ty cos wt 
Po 0 x=Ly Y ne 
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By eamspection, the solution for Pp, must be 


Py L, = k CoS kev! 
a et ieee Sioa oso 
Poo y x Se esc 
L. 7 k cos KL" x ony 
+ M mn > io Sin k "L cos oo cos wt (5.3) 
ye > aad y 
where kK! satisfies the equation 
| oe 
ky ae C54) 
and ae Satisfies 
rT 2 Tee ae Ww eZ 
Cee ee, & co oe (B.8) 


Solution for p, allows the creation of Po c Py + Spy 


~= 





Po ' Ly - k _ cos Lys 
5 = M cos KY cos wt -— 6M jae WE We ot Ste ey COS WE 
xc VY >,6 x xX 
0-0 
L, = KY cos k"x On 
+ 6M CL z jae a Cos ae cos wt ° (5.6) 
Y x x y 


~ 


Let the cavity dimensions be designed such that the (0,2,0;2) 
and the (1,1,0;2) modes will be degenerate. This requires 


chat 
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and 


= = 738 (5.7) 


for the desired degenerate modes. When (5.7) is inserted 


mreo (5.6), the result is 





0 1 er Komen COS meh 
= M cos ky cos wt - 6M — ie SLMS SEG Wt 
Pye v3. 7 
0 
Ke cos ki] 
te a cos ony COS SUE ° (5.8) 
eu ad x “x y 
This must be squared when substituted into (3.11). 
Recalling that 
(3.6) 


Po = Po + Spy ; 


26 





then the non-linear term has the form 


ies 2 Po P; ee il 
(5) = (—Sy) + 26(—— Sy) (5) ++ 8 (>) ise 
PAC pc oye: Dua O.c : 
oO 0 0 0 0 0 0 0-0 


The third term on the RHS above is of order Se and may 


therefore be neglected. The first term becomes 


0) o6=M cos’ ky SOs" DieLe 





Using the trigonometric identity 
cone a = > [fl + cos 2a] 


produces the following for Py 


Dowie ic 
= [1 + cos ae pie reose Zu) t | 


Since the only significant contribution will come from 
terms near resonance, those terms not involving a product 


of cos eK Y cos 2wt may be neglected. Thus 





Bonet Sah 
( 5) = 7 cos 2k y eee Ae we (CS. 109) 
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Similarly, keeping just those terms which will lead to 


resonance produces 


Pp Pp a k 
26 zg a =- See fly) a ee, Sa cos k_ycos2wt 
P,c P,c e Oe) aim teh Y 
0 0 0-0 x X 
k COstier x 
+- 

+: a SNe (2u)* — $Y, —— cos k_y cos 2uwt 

V3 xX Sin KL, 


(ie le) 


Micememncgmtme results of (5.9 —- 5.11) into (3.11) gives 


P 
Dene + Be: 
Co Be —2 =o Mw cos ky cos 2wt 
0-0 


k cos k_'x 


= + 6 Ac Din) 7 = a ee Os Vy CeCOS eit 
¥3 Ie oe lao eae 16 Y 
x x "xX 
k cos k._' x 
5 x 6 M (nie = >. = cos Ky CostZut. 
i ee ee stn kL, Y 


(5212) 


The terms cos Kk 'x and cos kx must be expressed as a cosine 
series of the form given by (4.7). When this is done, and 
only those components of argument Tx/L,, are retained, then 
all terms driving the (1,1,0;2) mode will be isolated, as 


desired. The results are 


ko! = w/L (5. 23a 


20 





TT 





cos kL 'x = 0.535 cos mee (5.135) 
x 
sin Kk, Ly = 0.971 , (Ses) 
KS iv3 7 (5.13d) 
L 
y 
cos k,"x = -3.677 cos — ; (5.13e) 
Xx 
and 
Sialicl kL, =i. (5) oe) 


where use was made of the identities 


Sinh A 


il 
}- 


Sie aN 


and 


cank A 


il 
Nol e 
OD 
| 
4) 


When the results of (5.13) are inserted into (5.12), the 


result is 


Pied ey ee 
Cy Nn 72 = M~ w cos aS cos 2wt 
Ome) 
- §61.83 a ie cos =~ cos _ cos 2wt, (5.14) 
a vy 
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where the first term on the RHS is called the finite 
amplitude term, and the second term on the RHS shows the 
effects of the boundary irregularity on the first order 


non-linear perturbation correction. 


Bees LHE SOLUTION FOR Py 

Appendix I shows the general solution for a function 
operated on by c ¢ ie Using those results, the solution 
icone Py for this example becomes immediate: for the finite 
amplitude term, designated by subscript f, 


an  - vy+1 Qe See Be 


re (Dye 


2 


( M ae cos BR cos (2wtt+B_) , (5.15) 


alee 


and for the boundary irregularity induced term, designated 


DyaeeSOOScript p, 


P Csi 
(——>) = -1.83 6 ae Me Wy? 5 cosk  xcosk ycos (2wttB.,) . 
Pog P (2w) Y 
(Se1:6)) 
When (5.15) and (5.16) are combined, they form P the first 


pe 
order non-linear perturbation correction to the pressure. 
Py 1s a non-linearly generated term which results from 


5 


the existance of a 010 standing wave in the cavity. Py 


Pp 
exists because the irregular boundary forces a non-linear 


interaction between the correction term Py and the 010 mode. 
It should be noted that there exists an additional 


correction term to the total pressure. This is the classical 
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correction to the 020 standing wave which exists ina 

cavity with an irregular wall. This additional correction 
term will be on the same order of magnitude as Py wae 

and must be taken into account when finding the total effect 
of the irregular boundary on the pressure. The following 
discussion cf results will consider only the correction term 


generated by this example. 


Cay RESOULTS 
In order to best illustrate the magnitude of the effect 
of the non-familial mode in this example, recall that from 


Or bS)eand (5.16) , 


P QO. Sin 8 
wk: = ye t+ ne a cos 2k y cos (2wttB-) 
0c (2w) Y 
0 0 
Om sin bs 
- 1.836 + Pb ___ PB Wc me cos kx cos k_y cos (2wttg_). 
(2) - © 


Dividing both sides of the equation by 


tee 


Z 
2 4 


allowing the quality factors Q- and 2, to be equal, and 


assuming measurements are made at one of the corners so that 


os k_ x = os k = os 2k = dl, 
cos ky Cc yy E yy 


S),1E 





gives 


P,/p e ‘ 

5 = sin Be cos (2wt+B -) =PieO30SIN eeeecos (Z2ZuUttsh). 
ytl M P P 
a od 


Let the LHS be designated by C. Additionally, introduce 


Q- 2, vt =F) 
Q¢ + O, W) 


H = 2 
aS a measure of the overlap of the curves obtained by 


plotting the magnitude of C versus a function of frequency 


for two cases: 


1 - 6 = 0, corresponding to no non-familial mode 
existance, and 


2 - 6 #0, corresponding to a degeneracy between the 
(Ot, 0-2 amd mle ilasOie2) modes). 


These curves, Figs. 2-6, allow a comparison of the magnitude 


of the P, term with and without the presence of a non- 


ul 


familial mode. The function of frequency on the abscissa 


is 

F = cot Bes 
where 8 is as defined by (I.1). Plots are provided for five 
different values of H. Additionally, for the case H = -l 
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there is a plot labelled H_ = -1 which shows the result of 
selecting a different corner for measurements, one 180° 
out of phase with the plot labelled as -~l1. 36 was arbi- 
Eiaraly chosen as 0.1. 

From the graphs it is clear that the presence of a non- 
familial mode does affect the total solution for the pressure 


by an amount proportional to 6, the magnitude of the boundary 


perturbation. 
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Viz CONecEVUSIONsS 


Non-linear theory has been applied to standing waves 
in a rigid-walled rectangular cavity with a perturbed 
boundary in an attempt to find one possible mechanism for 
the excitation of resonant modes other than those belonging 
to the family of the driven mode. It was observed that such 
an excitation does exist if the boundary perturbation and 
the dimensions of the cavity are favorably chosen. If the 
fractional perturbation of the boundary is given by 46, then 


the correction terms are of that same order. 


Sie, 





Ae een 


It is desirable to show that when ea 33 operates 
Sma FuNCtIOn of the form cos kx cos KY cos k,2 COS (UG )ae, 


the result will be of the form 


2 
ee 
Q 





= sane cos kx cos Kol cos k 2 cos (Htto—B)) 


where 8 is an angle defined such that 


i _ 4p ee 
5 (ae (=) esr BS - (ciela} 


Begin by recalling that, by definition, 


20 ig ae lle 
eer = oS." eae se (ee) 


Porrorming this operation on a function of the form 

cos kx cos KY cos k,z cos(wtt$) produces three terms, all 
with the common factors cos kK x cos KY cos k 2 . Canceling 
these like terms from both sides of (1.2), and using the 


fact that 


requires that 
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fey k cos (wt+¢) + Wcos (wt+d) + Sin(Wt+¢?) = 


2 
ron <= cos (Wt+o-6) 


Om, collection terms, 


C 2,2 2 


2 
2 p W W i 
w {l - ] cos (Wt+$d) + — Sin(wttd) = —~-— cos (wt+¢o-8) 
2 @) > QO sin B > 


Rearrangement of (I.1) gives 


= cos 86 
@ “Sie sa” 


and substitution of this result into (I.3) requires 


2 2 
cos (wt+¢) + oF Sin(wt+¢) = at 


iL 
Sin 8B 





cos (wt+¢-8) 





Canceling the common factors Haye , and multiplying both 


Sides by sin g then requires that 


cos 8cos(wttd) + sin 8 Sin(wttd) = cos(wtt+d-B) . (I. 4) 
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Using the trigonometric identity 
Gosh cos 5B + Sin A sin Bb) —~ ces(A = 8B) 


immediately shows (I.4) to be satisfied. Thus a general 
Porucron to ct when operating on a function of a 


specific form has been shown: 


a cos k,x cos Ae cos Kz cos (wtt+o>) = 


Z 





W il: 
ea—_—= +O— 
7 Same cos kx cos Magi cos k 2 cos (wtt+¢o-8) 
CA 
where tan B = 5 fit: ( 3 aa 
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